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Abstract 

We reconsider the bounds on non-standard neutrino interactions with matter which can be de- 
rived by constraining the four-charged-lepton operators induced at the loop level. We find that 
these bounds are model dependent. Naturalness arguments can lead to much stronger constraints 
than those presented in previous studies, while no completely model-independent bounds can be 
derived. We will illustrate how large loop-contributions to four-charged-lepton operators are in- 
duced within a particular model that realizes gauge invariant non-standard interactions and discuss 
conditions to avoid these bounds. These considerations mainly affect the 0(10^^) constraint on the 
non-standard coupling strength Se^, which is lost. The only model-independent constraints that 
can be derived are 0(10"^). However, significant cancellations are required in order to saturate 
this bound. 
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I. INTRODUCTION 



Neutrino non-standard interactions (NSI) were originally proposed |l| as a mechanism 
to produce neutrino flavour conversion in matter and considered as a possible explanation 
for the solar and atmospheric neutrino deficits . Although it is now 

clear that such NSI cannot fully account for the observed neutrino fiavour conversion, the 
increasing sensitivity of neutrino oscillation experiments to sub-leading effects has triggered 



a new interest in them anc 
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'erence with neutrino oscillations at present {e.g., K2K, 
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VjT) and future (e.g., SuperBeams, /?Beams or 
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301 1) facilities. In particular. 



the determination of the leptonic mixing angle 6*13 could be severely affected by degeneracies 



with the non-standard parameters 
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33|. 



Neutrino NSI can be described by effective four-fermion operators of the form 



(1) 



where / and /' are charged fermions with the same quantum numbers and L, R represent 
the chirality. For the rest of the paper we will denote with u the left-handed neutrinos, £ 
the left-handed charged leptons, L the left-handed lepton doublets and E the right-handed 
charged leptons. For NSI of neutrinos with normal matter, f = f can be either an electron, 
an up-quark or a down-quark. While strong experimental bounds are present on the cor- 



responding four-charged-fermion interactions, (in'y^iFi)( .fr,,Rni,..fi p), neutrino NSI are much 



less constrained. Constraints on them have been derived in Ref. |34l . l35l . l36l | making use of 
both tree level and one-loop processes. Here we reconsider the one-loop bounds, focusing on 
the necessity of implementing the neutrino NSI in a gauge invariant way in order to obtain 
a gauge independent result. 

The simple promotion of the neutrino fields in Eq. ([T]) to lepton doublets in order to 
construct a gauge invariant operator would imply that the strong bounds stemming from 
fiavour violating four-charged-fermion processes will apply to neutrino interactions as well. 
In order to avoid this and allow large NSI, the simultaneous presence of^tree level fiavour 
violating four-charged-fermion interactions must be forbidden 37], |38 



39( 1 . However, even 



if this requirement is satisfied, these interactions can be generated at one loop from the 
operator of Eq. ([T]) via a W exchange between the neutrino legs. In Ref. 34] this has been 
exploited to set bounds on some e. Notably, an 0(10""^) bound on eljj^ was derived through 



FIG. 1: One- loop contribution to the four-charged- fermion vertex arising from the operator of 
Eq. dl]) via W exchange. 



loop contributions to the decay jj, 3e and the /i — e conversion in nuclei. Consequently, the 



non-standard coupling strength elf was ne 
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ected in a large number of studies (see for exam- 
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44l|). However, as we will show, when gauge 



invariance is imposed, only naturalness arguments can be invoked. These arguments can 
lead to even stronger constraints on the NSI strength, but no completely model-independent 
bounds can be derived. The relevant diagram contributing to /i — > 3e is depicted in Fig. [H 
The computation of this diagram using only the operator of Eq. ([1]) with / = /' = e renders 
a gauge dependent result due to the fact that the operator itself is not gauge invariant. 
Thus, the necessity of considering a gauge invariant formulation of NSI is manifest. In the 
following, we will consider gauge invariant realisations of NSI both when the operator of 
Eq. ([T]) is realised from a dimension-six {d = 6) and from a d = 8 operator. We will discuss 
here the NSI with charged leptons, i.e., f = f = i. However, similar arguments as those 
presented here are applicable in the case of neutrino NSI with quarks. While gauge invari- 
ance has to be carefully taken care of in the cases mentioned above, one-loop processes like 
the W and Z decays (such as the ones depicted in Figs. 5-7 of Ref. [3^) do not contain 
gauge boson propagators in the loops and the bounds derived from them still apply. 



II. d = 6 REALISATIONS 



A. Flavour Antisymmetric Operator 



There is only one gauge invariant d = 6 effective operator that can induce the effective 
interaction of Eq. ([1]) in a direct way while avoiding the generation of four-charged-fermion 



operators. This operator is of the form: 



Ol = {^^T2K){Ll^TM . (2) 

where L'^ = CL^ and C is the charge conjugation operator. It can be generated upon 
integrating out a heavy charged scalar sing let with lepton number violating couplings to the 



lepton doublets, as in Ref . 45|, |46|, |47|, |48|| . This operator induces NSI only for leptons and 



with a very characteristic flavour structure^: 
This structure implies the antisymmetry relations 

0/3 _ 7/3 _ aS _ 7<5 ( A\ 

Thus, for each process of the type £ 3£' (i.e., a decay of a heavy lepton into three lighter 
leptons), there will be four contributing diagrams at one loop. For example, the process 
fi^fi^e^ will receive contributions from e^'^, e^f^ and e^^. In the approximation 
of massless fermions, using the antisymmetry relations of Eq. (jll), the four contributions 
exactly cancel. Thus, any loop contribution from d = Q operators to the four-charged- 
fermion interactions will be suppressed at least by a factor 0{mj/M^), where me is the 
mass of the heaviest lepton involved in the process. As a consequence, the bound on the 
corresponding e will be increased by the inverse of this factor. However, the antisymmetry 
relations of Eq. (jl]) also imply that NSI with electrons that change the neutrino flavour are 
related to charged lepton flavour changing interactions, and thus stringent constraints can 
be derived in this case (see 0, Q). Note that /i — > 3e is forbidden for (i = 6 operators, 
since there are no e with the given symmetries which can contribute to this process. 



B. NSI from Non-Unitarity 

The second possibility to generate neutrino NSI avoiding four-charged-fermion interac- 
tions is in an indirect way via the dimension six operator 

Of'^ = -{Lo^H)ip{H^Lp), (5) 

^ NSI with right-handed fields avoiding four-charged-fermion interactions can be reahsed through higher- 
dimensional operators as we will discuss in the next section. 
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where H is the Higgs doublet (we choose the hypercharge of H to be 1/2) and H = iT2H*. 
This operator induces non-canonical neutrino kinetic terms. After diagonalising and nor- 
malising them, a non-unitary leptonic mixing matrix is produced and, upon integrating out 
the W and Z bosons, neutrino NSI are induced. The tree level generation of this operator, 
avoiding similar contributions to charged leptons (see, e.g., Ref. |50l]), involves the addition 
of Standard Model-singlet fermions (right-handed neutrinos) which couple to the Hi ggs and 



lepton doublets via Yukawa couplings like in the standard seesaw model 
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This second realisation of NSI at c? = 6 is a 



a non-unitary mixing matrix (see Refs. 38 
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so quite constrained due to the effects of 



561]). The bounds on the SepL element are 



particularly strong due to the enhancement of the fi —>■ e'-f process if the GIM mechanism is 



not realised given the non-unitarity of the mixing matrix. In Ref. [38j, a bound of |ee/x| < 
5.9 X 10-^|ra„/rae - 1| for Mjy^ > Mw or \eef,\ < 9.1 x 10-^|ra„/ne - 1| for Mjy^ < Mw 
was computed, where n„ (ng) is the neutron (electron) density in matter. Notice that, since 
n„ ~ He, this means an additional suppression of (9(10^^). These bounds are already much 
stronger than the loop bounds discussed here and thus we will not consider this possibility 
further. 

III. d = 8 REALISATIONS 

Dimension eight realisations of the NSI offer more freedom. Gauge invariant d = 8 
operators can be generated by adding two Higgs doublets to the four-lepton operators. The 
vev of the Higgs field can then be exploited to break the SU{2) symmetry between the 
charged leptons and the neutrinos and can thus induce neutrino NSI avoiding their charged- 
lepton counterparts. A basis for operators involving two Higgs doublets, two left-handed 



lepton doublets anc 



is given in Refs. 



37 



;wo fermions, either left or right-handed (for the matter components). 



39| 



Oll^H = {Lf,YL^){Lsi,L,){H^H), (6) 

OTlh = {LpYTL^){Lsi,rL,){H^H), (7) 

Olfn = {LpYLa){LslpTL,){HWH), (8) 

Olll = {LpYrK){UlpL,){H^fH), (9) 

Of^H = i-ie'^"'){Lf,Yr^L^){Lsi,T'L,){H^T^H), (10) 

0™ = {LpYL^){Es^,E,) [H^H] , (11) 



01]!h = {LpYrL^m^pE,) {H^tH) . (12) 

In addition to the two left-handed lepton doublets, the two last operators contain two right- 
handed charged leptons and the first five two additional left-handed lepton doublets. The 
generalisation to operators involving interactions with quarks is straightforward replacing 
these fields by their quark counterparts. Generically, after electro-weak symmetry break- 
ing (EWSB), these operators generate both neutrino NSI and non-standard four-charged- 
fermion interactions at tree level. In order to avoid the latter, the following conditions have 



to be met [agj^: 



^111 _ _/^313 inlll ^331 , ^133 , ^313 _ r. ^333 ('<o\ 

The second condition involves four operators of the basis and we can thus choose three 
independent combinations satisfying it in order to form a basis for the d = % operators that 
induce neutrino NSI avoiding four-charged-lepton interactions: 

^ = liOlYn - OH'h) = mT,L,){UT,Ll){H^H) , (14) 

B = {Oil], - OITh) = 2{Lp^,L^){LsH)Y{H^L,) , (15) 

C = [Of^u - Oflu) = 2{L^H)Y{H^L^){LslpL,) , (16) 

^ = Olll, = 2{Lf,^,L^){LsH)Y{H^L^) - 2{LpH)Y{H^L^){LslpL^) , (17) 

^ = {Ol'iH - ODin) = 2{L^H)Y{H^L^){Esl,E,) . (18) 

All of these operators have to be divided by M*^, where M is the new physics scale from 
which the NSI originate, and multiplied by a flavour- dependent coefficient: a°'^'^^ multiplies 
A, h°'^'^^ multiplies B and so on. Note that, apart from which is independent since it is 
the only one containing right-handed fields, the other operators are all related. Indeed C is 
obtained from B simply by rearranging the fiavour indexes, T) is equivalent to i3 — C after a 
Fierz transformation and again a rearrangement of two indexes and also A can be written 
in terms of combination of B operators with different fiavour indexes. Thus, the only new 
structure that (i = 8 operators offer to select neutrino NSI avoiding four-charged-fermion 
operators is the combination {LH)'^'^{WL) present in Eqs. f|T5|) - f|T8|) . 

We will therefore compute the loop of Fig.[T]for this new structure. From the combination 
{L'ypL){LH)Y{H'^L), "opening" the operators in components, it is easy to check which 



^ Notice that our relations differ from those in Ref. 
charge. 



39| due to different conventions for the Higgs hyper- 
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FIG. 2: The effective interactions induced by B. 

operators are generated together with the one in Eq. ([1]). Taking all of them into account, a 
gauge invariant computation can be performed. In this case, the generated e are completely 
independent and the effective interactions shown in Fig. [2] are all generated with the same 
strength. Out of these diagrams, diagram (a) will result in the effective neutrino NSI when 
the Higgs acquires a vev. There are now several one-loop contributions to the four charged- 
lepton vertex. In diagram (a), we can connect the neutrino lines with either a ly or a 0+, 
resulting in a conversion of the neutrinos into charged leptons. For diagrams (c) and (d), 
we can connect the charged Goldstones to the neutrino lines, and finally, for diagram (b) 
we can close the 0+-loop. It is now possible to check the gauge independence of the result 
explicitly. We have done this by performing the computation in the i?^ gauge and splitting 
the W propagator in the unitary gauge part and the ^-dependent part and checked that the 
gauge dependence introduced by the W and all the diagrams with the Goldstones cancels, 
as it should. The remaining gauge independent contribution is then the one from the W 
exchange with the propagator in the unitary gauge. 

The most widespread way of estimating bounds from loop processes of an unknown high 
ener gy t heory through its effective description is exploiting the logarithmic divergence, as in 
Ref. [34] • Indeed the coefficient of this term and that of the logarithmic running of the full 
theory are the same and the mild scale dependence is just an 0{1) correction. On the other 
hand, the finite and quadratic contributions of the effective theory are less reliable, since 
they depend on the matching with the unknown full theory. Moreover, these contributions 
can be fine-tuned away, while for the logarithmic running this can only be true at a given 
scale^. 

We find that, neglecting lepton masses, the logarithmic divergences coming from the two 

^ Unless the cancellation is implemented through an operator with exactly the same running. 
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parts of the ly-propagator exactly cancel at one loop.^ Then, also in the case of this d = 8 
operator, only very weak bounds on e can be derived through the logarithmic divergence 
using the decay of a heavy lepton into three lighter leptons. It should be noted that, even 
if the logarithmic divergence is not present, a quadratic divergence is. We will devote the 
next subsection to the physical interpretation of this quadratic divergence. 

We argue that the discussion presented here can be also applied to higher-dimensional 
operators where only Higgs doublets are added. In order to preserve gauge invariance in 
the NSI, it is necessary to include the effects of the internal Goldstone loops. In order to 
avoid specifying the underlying theory, we must therefore compute the loop diagrams of the 
effective theory in the unitary gauge, where the Goldstone propagators vanish. 

A. The quadratic divergence 

The computation of the loop with a W exchange between the neutrino legs of diagram 
(a) of Fig. [2] in the unitary gauge turns out to give a vanishing logarithmic contribution, 
while a quadratic divergence is present. It is easy to check that a similar diagram with a Z 
exchange plus a diagram where the physical Higgses of diagram (a) of Fig. [2] are closed in a 
loop give exactly the same quadratic divergence to the neutrino NSI operator. 

It is simpler to understand the origin of this quadratic contribution to both operators 
from the decomposition into singlet minus triplet of Eqs. ( fT5l) - (fT8l) : 

{LH)Y{H^L) = ^ [{LYL) {H^H) - {LYtL) (HWH)] . (19) 

If the Higgs legs are closed in a loop, the resulting tadpole will contribute with a quadratic 
divergence only through the singlet term. The corresponding divergence is not present for 
the triplet term since it involves a trace of the r matrices, which vanishes. This means 
that the quadratic divergence will be proportional to the singlet combination (L^^L) and 
thus will be equal for the neutrino NSI and the four-charged-fermion interaction. Since the 
logarithmic contributions to the process vanish, we can now try to estimate the constraints 
that this quadratic divergence implies. We will consider these quadratic contributions for 
the operator of Eq. f|T8l) involving the right-handed leptons E, but similar discussions also 
apply to Eqs. (fT5|) - (fT7j) . as well as to interactions with quarks, while the contributions cancel 

^ Nevertheless, they might be present at the two loop level. 
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for Eq. (fT^ . However, as we will show later, for the left-handed fields antisymmetries similar 
to the ones discussed at ci = 6 can be considered in order to avoid the bounds. 
At tree level, after EWSB, the operator E generates neutrino NSI with strength 

On the other hand, if we parametrise the four-charged-lepton interactions as 

£cF = -1^Gpe^^^l{U^/^){W,^^E,) , (21) 

then B^Qp]^ = 0. However, the contribution from the quadratic divergence arising at one 
loop is equal for the neutrino and the four-charged-lepton NSI. Adding this, we have: 



^/3a — o ;i^2 ( ;i^2 + a^2o_5 ) (2^) 



2 M2 M2 87r2_ 

^CF,R - o /\^2/\#2o„2' V^"^/) 



2 M2 M2 87r2 

where = 0{1) depends on the UV completion of the full theory, and A < M is the 
scale at which new physics appears and the effective theory is no longer valid. Thus, it is 
clear that no model independent bounds can be derived from this quadratic contribution, 
since assumptions on the sizes of A and k have to be made. For example, the high energy 
completion of the theory could be such that k = due to some significant fine-tunings and 
then no bounds would stem from this process. On the other hand, naturalness arguments 
can be invoked to argue that, in absence of significant fine-tunings, k = 0{1) and A should 
be at least as high as the electroweak scale. With these assumptions, bounds of the same 



order as the ones derived in Ref. SJ] will be recovered. On the other hand, if no new physics 
is present between v and M, A could be identified with M and very stringent constraints 
could be derived. 

Below we will consider, as an example, the computation of the quadratic contribution 
in a complete theory whose low energy effects are described precisely by the operator of 
Eq. f|T5|) . without the need of fine-tuning in order to cancel similar operators contributing 
to four-charged-fermion processes at tree level. In this example, A = M and k = 1/2. 
Thus, even if only neutrino NSI are induced at tree level, the loop contribution only has a 
suppression of Stt^ ~ (9(100), which could dominate the 0{v^/M'^) tree level contribution 
(unless M < 1 TeV) and the four-charged-fermion operator would be induced with a strength 
similar to that of the neutrino NSI. This would imply strong bounds on the neutrino NSI 
and the main motivation for considering d = 8 operators would be lost. 



We will now discuss how some antisymmetries between the e parameters could avoid 
these corrections, in a way similar to the d = Q antisymmetric realisation. As for = 6, this 
is only possible for the case in which the matter fermions are left-handed leptons. To study 
this, it is more convenient to change the operator basis to a basis where these symmetries 
are manifest: 

- - (pa'^Mfiip^^)] , (24) 

+ |0+n[(z7^7^z/,)(47p^^) + (z>57"z/^)(Vp4) 



2 



X 



+2 ^U^{pf,^Pv^){yapV^) + + • • • . (25) 

-S + C _ 
2 ~ 

(10+1' - \U^W^in:){v,^pV,) - (47%)(^/37p^a)] + • • • , (26) 



where . . . represents terms proportional to 0o'/'+ or 0o0+) which will not contribute to our 
one-loop computations as long as we consider the leptons to be massless. 

In the Feynman gauge, the quadratic divergences are completely determined by the loops 
of the Goldstones and Higgs fields. Indeed, since the quadratic divergence is independent of 
the mass propagating in the loop, the terms |0op (via the Higgs and neutral Goldstone) and 
(via the charged Goldstones) will give the same contribution. It is thus evident, by 
looking at Eqs. fl24l) - fl26l) that the operators S and A will contain the quadratic divergence, 
while X will not, since the two contributions cancel. 

Also, a generic coupling 

^a/37<5 can be decomposed as 

where 

^a/37(5 _ ^iPaS _ ^aS'yfS _ ^-ySafS _ ^(^0/375 _|_ ^7/3a(5 _|_ ^«57/3 _|_ ^-ySafS^ {'28) 
_ _^7/3ck5 _ _^a5-yl3 _ ^jSaP _ }i^^a(Sj5 _ ^-y/SaS _ ^aS-yfi _|_ ^7<5o/3^ (gg^ 

We see that s has the same symmetry as S, a has the symmetry of A, and x has the 
symmetry of X, and are actually the coefficients of S, A, and X. 
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FIG. 3: An example of a full theory inducing neutrino NSI, but not four-charged-fermion inter- 
actions, at tree level. The NSI are generated through the exchange of right-handed neutrinos Nr 
and a scalar doublet S. 



To illustrate how operators with different symmetries contribute in a different way to the 
quadratic divergence, we will take the theory from Section 5.1 of Ref. {5] as an example (see 
Fig. [3]). In this theory, the NSI are realised at tree level by two Higgs doublets selecting the 
neutrinos from two lepton doublets by taking a vev and the exchange of two right-handed 
neutrinos and a scalar doublet. For simplicity, we will here assume that the scalar has 
the same mass M as the right-handed neutrinos. The effective tree level operator is now 
essentially given by 

'^^{LpH)Y{H^L^){LsipL,) = ^(5 + A + X) (31) 

and thus contains all of the above mentioned operators, of which two have the quadratic 
divergence. At tree level, after EWSB, it generates neutrino NSI and four-neutrino in- 
teractions^ with strength = Notice that the introduction of right-handed 
neutrinos already implies that NSI will be induced at d = 6 through the deviations from 



unitary mixing. Consequently, the constraints derived from non-unitarity in Ref. [38|] would 
apply and the loop bounds would not be so relevant. However, we find this toy example 
useful to connect the quadratic divergence to a full theory in which it can be computed and 
matched in order to clarify its interpretation. Indeed, since this Standard Model extension is 
renormalizable and the corresponding four-charged-fermion operator does not appear at tree 
level, all diagrams are actually finite and we will be able to calculate them unambiguously. 
When computing the loops of the Higgs and Goldstones in the full theory, we will assume 
that they are essentially massless compared to the heavy mass scale M. We will also assume 



These are generally hard to constrain directly and their effects are usually relatively weak, see, e.g., 
Refs. [silsilsSl. 



11 



that the external momenta are neghgible. With these approximations, the loop contribution 
is 

+ iiyaYl^fB)ii'ylp^5) + iiy^Yl^5)iL7p^/3)] ■ (32) 

It is easy to check that this is proportional to 5 + as anticipated. We note that this 
operator can be obtained from Eq. fl3Tl) simply by replacing \(f)o\'^ and by the factor of 
M^/ (327r^) coming from the loop integral. Thus, the complete one-loop four-fermion vertexes 
are given by replacing |0op with f^/2, in order to determine the tree level contribution, and 
replacing both |0oP and by M^/(327r^) in order to determine the loop contribution. 

From Eq. fl52]) we see that the following interactions are generated, all with similar strength: 
four-neutrino interactions, neutrino NSI, and four- char ged-fermion interactions. However, 
since the four-charged-fermion and the four-neutrino interactions are completely symmetric 
under the exchange of flavour indexes, while neutrino NSI are not, the remaining terms for 
neutrino NSI and four-charged-lepton interactions are: 



e 



«/3 



2 r „,2 1 



(33) 



Here, e^'F^f defined through 

£cF = -2V2GFe''c^fJ{laY^p){i,lp(-6) (35) 

Thus, if Svr^f ^ <^ M^, both the neutrino NSI and the four-charged-lepton operators will be 
dominated by loop effects. 

As we have just seen in the above example, a physical meaning can be attributed to the 
quadratic divergences obtained in the effective theory. Essentially, they can be regulated 
by reinserting the missing propagators of the heavy particles in the full theory, leaving a 
contribution to the effective NSI which is suppressed by f^/(87r^M^) instead of the tree 
level V* /M^. In addition, the completely symmetric contribution from s^^^'^ will generate 
an additional four-charged-lepton operator at the one-loop level. Thus, the s^'^^'^ of a model 
could be severely constrained by the strong bounds on decays such as /i ^ 3e, while the 
parameters a^^^"' and x^^^"' only contribute to the neutrino NSI and cannot be constrained 
from these processes. However, it is challenging to build a full theory which generates the 
antisymmetric couplings, but not the symmetric one, in a natural way. 

12 



An important caveat: in the above example of a full theory, we assumed that the masses of 
the different heavy particles were the same. In general, the tree and loop level contributions 
to the es may be different functions of the mass ratios and couplings, meaning that it could 
be possible to fine-tune these functions in such a way that the loop-contribution is zero, 
while still maintaining a non-zero tree level contribution. 



IV. SUMMARY AND CONCLUSIONS 



We have reconsidered the bounds on neutrino NSI from one-loop processes. We have 
shown that, in order to have non-ambiguous bounds, a gauge- invariant realisation of the 
NSI must be considered. We explicitly studied d = 6 and d = 8 operators and have shown 
that, in both cases, the logarithmic divergences of the one- loop contributions are suppressed 
by the factor mj/M^, which severely weakens the bounds. 

In particular, for (i = 6 operators, the anti-symmetry relations that arise as a consequence 
of the requirement of the absence of four-charged-fermion interactions at tree level force the 
one-loop processes to be zero in the absence of leptons masses. 

In the d = 8 case, we have shown that the loop processes involving NSI should be 
calculated in the unitary gauge in order to obtain a gauge invariant result if only the W 
exchange diagram is considered. In this way, the Goldstone loops present in gauge invariant 
realisations of the NSI with extra Higgs doublets are automatically taken into account. The 
result is that the logarithmic divergence is proportional to the factor mj/2M^. However, 
a quadratic divergence is present and can be exploited to set bounds on NSI. The use of 
the quadratic divergence in such a way implies model-dependent naturalness assumptions, 
in particular on the coefficient of the divergence and the size of the cut-off scale. 

For a coefficient k ~ 0{1) and a cut-off A of the order of the electroweak scale, the bounds 



presented in Ref. 



3J] are recovered. Pushing the cut-off scale to M, where the effective 



operators are generated, effectively assuming that no new physics appears between the 
electroweak scale and M to cancel the quadratic contribution, implies that the loop processes 
can dominate over the tree level contributions inducing four-charged-fermion interactions 
of a strength similar to that of the neutrino NSI, unless M < 1 TeV. This allows the 
derivation of strong bounds on the latter, but only based on naturalness arguments and 
not model-independently. All these considerations apply to NSI of neutrinos with both 
leptons and quarks. However, for neutrino NSI with left-handed leptons, we have shown 
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that, decomposing the NSI in parts with different ffavour symmetries, only the symmetric 
one contributes to the four-charged-fermion process. Thus, if this part is not present in 
the full theory from which NSI are realised, the loop constraints can be avoided and large 
neutrino NSI are still viable if generated via the antisymmetric couplings a and x. It remains 
an open question whether there are natural models which can realise this. Otherwise an 
extra fine-tuning beyond the one required at tree level would be necessary to cancel the loop 
induced four-charged-fermion operators. 

On the other hand, these naturalness arguments can always be evaded if one allows fine- 



tuning of the theory and (given that Refs. 



38 



39| have shown that large NSI avoiding 



four-charged-fermion operators already require a significant amount of fine-tuning at tree 
level) invoking naturalness arguments at the loop level will not make the model more natural. 
Thus, we emphasise the fact that no model-independent bounds can be derived from the 
loop processes studied here, but that the only viable models of large NSI avoiding four- 
charged-fermion interactions require significant cancellations not only at tree level, but also 
at the one-loop level. 

For practical purposes, in order to realise NSI without ad hoc cancellations, the bounds 



derived in Ref. 38| have to be respected or the neutrino NSI will be of the same order than 
the four-charged fermion operators. On the other hand, in a model independent approach 
considering possible cancellations both at tree and loop levels, the strongest effect is the 
loosening by three orders of magnitude of the bound on ell, since the constraints in the 
other flavours are dominated by tree level considerations {34]. These tree level constraints 
of (9(10^^) still apply to e{j^, but the stringent radiative bounds of 0(10^^) from ft 3e 
(/ = e) and /i — e conversion in nuclei {f = q), is lost. Given the strength of this bound, 
Sen = has been assumed for simplicity in many phenomenological studies. Therefore, it 
could be of interest to consider larger values for this parameter in order to determine its 
impact on future neutrino oscillation experiments. However, we would like to stress that 
this kind of large neutrino NSI would require significant fine-tunings at both the tree and 
loop levels. 

Acknowledgments 

We would like to thank S. Antusch, S. Davidson, A. Donini, B. Gavela, A. Hoang, 
A. Ibarra, V. Mateu, M. Papucci, C. Pena-Garay, S. Pozzorini, G. Raffelt, and I. Stewart 



14 



for discussions and in particular N. Rius and A. Santamaria for very useful and interesting 
discussions and for comments on this manuscript. M. B. and E. F.-M. would also thank the 
people at IFIC and University of Valencia for their kind hospitality during their visit. 

This work was supported by the Swedish Research Council (Vetenskapsradet), contract 
no. 623-2007-8066 [M. B.]. 



L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

S. P. Mikheev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 
E. Roulet, Phys. Rev. D44, 935 (1991). 

M. M. Guzzo, A. Masiero, and S. T. Petcov, Phys. Lett. B260, 154 (1991). 
G. Brooijmans, (1998), |hep-ph/9808498 



M. C. Gonzalez-Garcia et al, Phys. Rev. Lett. 82, 3202 (1999), |hep^ph7980953l l 

S. Bergmann, M. M. Guzzo, P. C. de Holanda, P. I. Krastev, and H. Nunokawa, Phys. Rev. 



D62, 073001 (2000), hep-ph/0004049 



M. M. Guzzo, H. Nunokawa, P. C. de Holanda, and O. L. G. Peres, Phys. Rev. D64, 097301 
(2001), |hep-ph/0012089j 

M. Guzzo et al, Nucl. Phys. B629, 479 (2002), ,hep-ph/0112310] 



Y. Grossman, Phys. Lett. B 359 (1995) 141, hep-ph/9507344 



T. Ota and J. Sato, Phys. Lett. B 545 (2002) 367, hep-ph/0202145 



A. Friedland and C. Lunardini, Phys. Rev. D72, 053009 (2005), 'hep^ph/0506143 



N. Kitazawa, H. Sugiyama, and O. Yasuda, (2006), |hep-ph/0606013 



A. Friedland and C. Lunardini, Phys. Rev. D74, 033012 (2006), hep-ph/060610l| 



M. Blennow, T. Ohlsson, and J. Skrotzki, Phys. Lett. B660, 522 (2008), [hei>ph/0702059 
A. Esteban-Pretel, J. W. F. Valle and P. Huber, Phys. Lett. B 668 (2008) 197, 0803.1790. 
M. Blennow, D. Meloni, T. Ohlsson, F. Terranova and M. Westerberg, Eur. Phys. J. C 56, 
529 (2008), 0804.2744. 

M. C. Gonzalez-Garcia, Y. Grossman, A. Gusso and Y. Nir, Phys. Rev. D 64 (2001) 096006, 
hep-ph/0105T59l 



A. M. Gago, M. M. Guzzo, H. Nunokawa, W. J. C. Teves and R. Zukanovich Funchal, Phys. 
Rev. D 64 (2001) 073003, |hep-ph7ol05T96[ 



15 



[20] P. Huber and J. W. F. Valle, Phys. Lett. B 523 (2001) 151, |hep-ph/0108193 



[21] T. Ota, J. Sato and N. a. Yamashita, Phys. Rev. D 65 (2002) 093015, hep-ph/0112329 



[22] M. Campanelli and A. Romanino, Phys. Rev. D 66 (2002) 113001, |hep-ph/0207350 



[23] M. Blennow, T. Ohlsson and W. Winter, Eur. Phys. J. C 49 (2007) 1023, hep-ph/0508175 



[24] J. Kopp, M. Lindner and T. Ota, Phys. Rev. D 76 (2007) 013001, [hep-ph/0702269 

[25] J. Kopp, M. Lindner, T. Ota and J. Sato, Phys. Rev. D 77 (2008) 013007, 0708.0152. 

[26] N. C. Ribeiro, H. Minakata, H. Nunokawa, S. Uchinami and R. Zukanovich-Funchal, JHEP 

0712 (2007) 002, 0709.1980. 
[27] ISS Physics Working Group, A. Bandyopadhyay et al, (2007), 0710.4947. 
[28] N. C. Ribeiro, H. Nunokawa, T. Kajita, S. Nakayama, P. Ko and H. Minakata, Phys. Rev. D 

77 (2008) 073007, 0712.4314. 
[29] J. Kopp, T. Ota and W. Winter, Phys. Rev. D 78, 053007 (2008), 0804.2261. 
[30] M. MaUnsky, T. Ohlsson and H. Zhang, Phys. Rev. D 79, 011301 (2009), 0811.3346. 



[31] P. Huber, T. Schwetz and J. W. F. Vahe, Phys. Rev. Lett. 88 (2002) 101804, |hep-ph/0111224 



[32] P. Huber, T. Schwetz, and J. W. F. Valle, Phys. Rev. D66, 013006 (2002), [hep^ph/0 2020481. 
[33] T. Ohlsson and H. Zhang, Phys. Lett. B 671 (2009) 99, 0809.4835. 

[34] S. Davidson, C. Pena-Garay, N. Rius and A. Santamaria, JHEP 03, Oil (2003), 
hep-ph/0302093| 



[35] J. Barranco, O. G. Miranda, C. A. Moura, and J. W. F. Valle, Phys. Rev. D73, 113001 



(2006), |hep-ph/0512195 



[36] J. Barranco, O. G. Miranda, C. A. Moura, and J. W. F. Valle, Phys. Rev. D77, 093014 
(2008), 0711.0698. 



[37] Z. Berezhiani and A. Rossi, Phys. Lett. B535, 207 (2002), |hei>ph/0111137 
[38] S. Antusch, J. P. Baumann and E. Fernandez-Martinez, Nucl. Phys. B 810, 369 (2009), 
0807.1003. 

[39] M. B. Gavela, D. Hernandez, T. Ota and W. Winter, Phys. Rev. D 79, 013007 (2009), 
0809.3451. 

[40] A. Friedland, C. Lunardini, and C. Pena-Garay, Phys. Lett. B594, 347 (2004), 



hep-ph/0402266^ 



[41] A. Esteban-Pretel, R. Tomas, and J. W. F. Valle, Phys. Rev. D76, 053001 (2007), 0704.0032. 
[42] M. Blennow and T. Ohlsson, Phys. Rev. D 78, 093002 (2008), 0805.2301. 



16 



W. Winter, Phys. Lett. B 671, 77 (2009), 0808.3583. 
D. Meloni, T. Ohlsson, and H. Zhang, (2009), 0901.1784. 
A. Zee, Phys. Lett. B93, 389 (1980). 
A. Zee, Phys. Lett. B161, 141 (1985). 

M. S. Bilenky and A. Santamaria, Nucl. Phys. B420, 47 (1994), |hep-ph/9310302' 
R. Barbier et al, Phys. Kept. 420, 1 (2005), |hep-ph/0406039j 
F. Cuypers and S. Davidson, Eur. Phys. J. C2, 503 (1998), |hep-ph/9609487, 
A. Abada, C. Biggio, F. Bonnet, M. B. Gavela, and T. Hambye, JHEP 12, 061 (2007), 
0707.4058. 

P. Minkowski, Phys. Lett. B67, 421 (1977). 

R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

T. Yanagida, In Proceedings of the Workshop on the Baryon Number of the Universe and 

Unified Theories, Tsukuba, Japan, 13-14 Feb 1979. 

M. Geh-Mann, P. Ramond, and R. Slansky, Print-80-0576 (CERN). 

D. Tommasini, G. Barenboim, J. Bernabeu, and C. Jarlskog, Nucl. Phys. B444, 451 (1995), 
|hep-ph/9503228[ 

S. Antusch, C. Biggio, E. Fernandez-Martinez, M. B. Gavela, and J. Lopez-Pavon, JHEP 10, 
084 (2006), lhep^ph/0607020 



M. S. Bilenky and A. Santamaria, (1999), hep-ph/9908272 



M. Blennow, A. Mirizzi, and P. D. Serpico, Phys. Rev. D78, 113004 (2008), 0810.2297. 



17 



